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INTRODUCTION 


Sadlation forms a significant fraction of the heat transferred to a probe 

that enters the atmosphere of Jupiter or Saturn. In reality, the design of the 

heat shield for such a probe is controlled by the radiation heat transfer. At 

present, the United States does not have experimental facilities capable of 

producing the extreme thermodynamic conditions necessary to test the heat shield 

designs for Jupiter or Saturn entry. Thus, the heat shields must be designed 

based upon analytical models of the flow field and heat transfer phenomena. 

The theoretical and numerical models must be as accurate as possible. In order 

(1 2 ) 

to improve the heat transfer predictions, Leibovitz and Kuo * and Tiwari and 

(3 A) 

Szema * Investigated the influence of finite rate ionization in the inviscid 
part of the stagnation shock layer on radiative heating of probes during Jupiter 

and Saturn atmospheric entry. Both investigative teams considered cold wall non- 
ablating models for the probe shock layer. However, they reached opposite con- 
clusions as to how the ionization rate assumption affects the radiative transfer. 
(3 4) 

Tiwari and Szema * contended that the radiation heating at the probe stagna- 

(1 2 ) 

tion point increased, whereas Leibovitz and Kuo ’ predicted that it decreased 
relative to the results obtained when the ionization rate was assumed to be in 
equilibrium at the local thermodynamic conditions. 

Present State of Knowledge 

The first studies of finite ionization rate effects on the radiative heat- 
ing of probes entering the atmospheres of Jupiter, Saturn, and Uranius were pre- 
sented by Leibovitz and Kuo^^*^^ and Howe^ ^ . These studies Indicated that the 
assumption of a finite ionization rate predicts significant reductions in the 

radiative heating as compared to equilibrium ionization predictions. Leibovitz 
( 2 ) 

and Kuo found that a nonequilibrium ionization layer with reduced radiative 
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emission exists just behind the stagnation point shock wave during part of the 

entry into model atmospheres of Jupiter and Saturn. They predicted that the 

radiative heating rate, calculated using the finite ionization rate assumption, 

could be up to 40% lower than if the heating rate were calculated by assuming 

that the plasma ionization takes place in equilibrium. However, the reduction 

in total probe heating, which was calculated by integrating both the radiative 

and convective heating rates over the atmospheric entry time, was less than 15%. 

The low reduction in total heat transfer occurred, because, when the finite 

ionization rate radiative fluxes were very low (compared to the equilibrium 

ionization rate radiative fluxes), the heating was largely due to convection. 

( 2 ) 

Leibowitz and Kuo calculated the radiative flux to the probe for both 

the equilibrium and finite ionization rate assvimptions by using the results of 

( 6 ) 

the hydrogen-helium radiation computations of Stickford . They computed the 
instantaneous values for the radiative flux at the thermodynamic conditions 
obtained from both the finite and equilibrium ionization rate assumptions at 
each point along the trajectory and Integrated the results to obtain the total 
radiative heating. 

(3 4) 

Tiwari and Szema ’ also considered the effect of finite ionization rates 

on the radiative heating of probes entering the Jovian atmosphere. Their results 

showed that the assumption of a finite ionization rate increased the radiative 

heating rate by about 10% as compared to the heating rate obtained by assuming 

equilibrium ionization. They used a fundamental fluid mechanics model and the 

(1 2 ) 

same ionization rate constants that Leibowitz and Kuo * used; however, they 
employed a detailed 58-step, spectral, absorption coefficient model for the hy- 
drogen-helium plasma to calculate the radiation transport. 

Zoby and Moss^^^ conducted a preliminary investigation of the thermal en- 
vironment and heating rates for a probe entering the atmosphere of Saturn. They 
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used the numerical code developed by Tiwari and Szema to determine the effect 
of a finite ionization rate in the inviscid part of the stagnation shock layer 
on the heat transfer to the probe. Their results indicate that the assumption 
of a finite ionization rate would significantly increase the radiative heat 
transfer rate at the probe stagnation point as compared to equilibrium ioniza- 
tion and that the finite ionization rate effects on radiative heating would be- 
come negligible away from the stagnation point. However, they did not allow the 
electrons and atoms to have different temperatures. 

These previous investigations show that the assumption of a finite ioniza- 
tion rate influences the radiative transfer and, hence, the design of the heat 
shields of probes intended for atmospheric entry of major planets. The non- 
equilibrium ionization layer just behind the shock wave influences the radiation 
transfer at or near the stagnation point. However, the radiative heating for 
the assumption of a finite ionization rate has been shown either to increase 
or decrease relative to its value for equilibrium ionization for roughly the 
same physical conditions in two independent investigations. This contradiction 
is investigated in the present study. 

The present study is limited to an investigation of the underlying reasons 

(1 2 ) 

for the opposite conclusions reached in the work of Leibowitz and Kuo * and 

(3 A) 

Tlwarl and Szema * for hydrogen-helium shock waves with a cold nonblowing 
wall boundary condition at the probe heat shield. The effects of ablation 

and the ablation layer gases are not considered. The study la limited to 

the stagnation shock layer. Appendix A contains a survey of the present 
state of knowledge of nonequilibrium shock wave structure. 
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SHOCK WAVE STRUCTURE 


This section contains a short introduction to nonequilibrium shock wave 
structure. A detailed review is included in Appendix A. Figure 1 shows the 
shock wave structure for a typical atomic gas taken from reference 8. The 
figure is for a Mach 24 shock wave in argon; however, it is similar to 
a hydrogen-helium shock wave, because the hydrogen dissociates almost instantly 
behind the shock. The pressure in front of the argon shock is 1 cm Hg, and the 

temperature is 300 K. The parameters shown are all nondimensional with tempera- 
tures, Tj, » T^/ (16,828 K), i = a, e, (where "a” denotes atoms and "e" denotes 

electrons) and density p - p/(. 000229 gm/cm ) . The reference value of the degree 

of ionization is ■ 0.449. The parameter represents the ground state 
continuum optical thickness. The reference values of temperature, density, 
and oLq are the equilibrium values behind the shock wave. Note that the 
nondimensional parameters representing these variables all go to unity as the 
gas goes to equilibrium. The nondimensional atom— atom and electron— atom 
ionization rates are given by and respectively, where n^ " 

1.31 X 10^^ l/(cm^-sec). 

The gas behind the shock wave is in thermal nonequilibrium, with the heavy 
particles much hotter than the electrons (f^ > T^) . As one moves into the re- 
laxation region behind the viscous shock, T^ rapidly increases to about its equi- 
librium value and then slowly increases until it becomes equal to T^ and thermal 
equilibrium occurs. 

Initially, just behind the shock wave, the degree of ionization is increased 
by atom-atom collisions. When the electron population is sufficiently high, 
electron-atom collisions become the major producers of electrons. Each elec- 
tron produced via an electron-atom collision requires that the electron gas as 
a whole supply the necessary ionization energy. Thus, the electron gas is cooled 
by ionizing collisions with the atom gas. At the same time, the electron gas re- 
ceives energy by elastic collisions with the atom gas and thereby reduces T • 

Si 
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Throughout most of the relaxation region, the energy gained by the electron 
gas through elastic collisions with the atom gas Is almost balanced by the loss 
of energy through electron-atom Ionizing collisions. This keeps the electron 
temperature roughly constant throughout most of the relaxation region. 

The argon atom was modeled as If It had only two electronic states: a 
ground state and an excited state as shown schematically In Figure 2. The 
single excited state of the model represents all the excited states of the atom. 
Thus, the radiation model involves only one bound-bound process and two bound- 
free processes. This type of two-state model has been widely used In nximerical 
shock wave structure studies. It generally predicts results that agree well 
with experiment. 

Figure 1 shows that the ground-state continuum optical thickness, T , has 

S 

a value of about 8 at the end of the relaxation region. This Implies that the 
ground -state continuum radiation emitted near the shock wave will not penetrate 
to the end of the relaxation region, because most of the radiation energy Is ab- 
sorbed within two optical path lengths of its emission point. The excited 
state continuum optical thickness is too small to be shown on Figure 1. Thus, 
excited state continuum radiation emitted just behind the shock can easily pene- 
trate to the end of the relaxation region. However, the excited state source 
function is very weak just behind the shock, because there are not very many 
free electrons available to recombine to an excited atomic state, and in the 
process emit radiation. 
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Figure 2, Schematic of the Radiation Absorption Coefficient for a Two-Level Atomic Gas. 




NONEQUILIBRIUM MODELING 


Collision Processes 

Leibowitz^^^ developed a reaction scheme describing the important col- 
lisional processes in hydrogen-helium ionizing shock waves that were modeled 
after results obtained from studies of argon ionization. Reactions are in- 
cluded for dissociation of molecular hydrogen, excitation of the electronic 
states of atomic hydrogen and helium, and ionization of the hydrogen and 
helium by collisions with atoms and electrons. The reaction scheme and rate 
constants are listed in Table 1. 

In the reaction scheme, the atomic hydrogen and helium are assumed to have 
two bound states: a ground state and a first excited state as shown schematically 
in Figure 2. The first excited state is a representative state, which represents 
all the excited states. The two-state assumption is a good model when the ground 
and first excited states are separated by a large energy gap as compared to the 
energy difference between the first excited state and the ionization energy. 

For these types of gases, the excited states and electrons have nearly the same 
energy. Thus, it is also assumed that the excited state and electron populations 
are in equilibrium with each other at the electron temperature. 

The dissociation rate for molecular hydrogen was obtained from shock tube 
investlgations^^^ . Hydrogen initially dissociates by collisions with itself 
and helium and, subsequently, by collisions with atomic hydrogen, ions, and 
electrons. For the thermodynamic conditions of interest in outer planet atmos- 
pheric entry, dissociation is completed before ionization begins, so that the 
dissociation and ionization processes are uncoupled. 

The ionization of atomic hydrogen is initiated by atom-atom collisions, 
which produce atoms with electronically excited states. The excited atoms are 
ionized by additional collisions. The reaction scheme assumes that hydrogen 
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Table 1 


Reaction scheme and rate constants for hydrogen- helium 
chemical nonequilibrium conditions, (From Leibowitz, 
Ref. 1) 



Reactions 

3 -1 -1 

Rate constants in cm sec mole 

1. 

H + e h'*’ + 2e 

1/2 

kj^ = 2.27E13 Tg exp(-157,780/Tg) 

2. 

He + e J He"^ + 2e 

1/2 

k 2 = 1.33E13 Tg exp (285, 160 /Tg) 

3. 

"k 

H + e H + e 

1/2 

k 3 = 4.11E13 ' exp(-116,010/Tg) 

4. 

* 

He + e He + e 

k^ = 2.24E13 exp(-232,030/Tg) 

5. 

H + H t H* + H 

1/2 

k 5 = 6.20E10 T ' exp(-116,010/T) 

6. 

H + He J H* + He 

1/2 

kg = 4.89E10 T ' exp(-116,010/T) 

if 

H 2 + He H + H + He 

k^ = 4.33E18 [1 - exp(-15E8/T^)] 

exp(-52340/T) 

8. 

H 2 + H 2 J H + H + H 2 

kg = 2.5 k^ 

9f 

H2+h:^H + H + H 

kg = 14.0 k^ 

10. 

H 2 + h'*’ H + H + H^^ 

^10 “ ^9 

11. 

H2 + e^H + H + e 

II 


All temperatures are ®K. 


^Reaction rates k? and modified from the expressions of Leibowitz by 
Tiwari and Szema^^>^\ 
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in the first excited state is produced by collisions with hydrogen (reaction 5) 
and helium (reaction 6) . Excitation to the first excited state is the rate 
limiting step for ionization. The excited atoms are rapidly ionized by further 
collisions. 

As the result of atom-atom collisions, the population of free electrons 
increases as one moves deeper into the shock layer. When the number of elec- 
trons reaches a certain magnitude, the electron-atom collisions become more 
important than the atom-atom collisions. Electron-atom collisions produce 
hydrogen ions by the same type of two-step process as the atom-atom collisions 
do. The atom is first excited by reactions 3 and 4 and then it is ionized. 

In summary, the reaction scheme developed by Leibowitz involves a two-step 
ionization process in which the excitation to the first excited state is the 
rate controlling reaction. In other words, the rate of excitation is finite 
in the model and the first excited state is not populated to its Boltzmann 
value. Thus, the Leibowitz model really allows for both nonequilibrium ex- 
cited state and electron populations. The finite rate ionization model also 

includes the assumption that the excited state population is in equilibrium 

(9) 

with the electron population at the electron temperature . Consequently, 
the model predicts low excited state populations just behind the shock wave, 
as compared to equilibrium chemistry solutions. 

Radiation Processes 

The general equations for the radiation absorption coefficient and source 
function are developed for a nonequilibrium atomic gas in Appendix B. The 
radiation absorption coefficient as given by Eq. 36, Appendix B, can be re- 
written for a two-level atom as follows; 

k^ = ng(l)aj^j[l - Zj^^exp(-hv/kT^)] + - Z2^exp(-hv/kT^) ] 

+ n^(l)aj^2t^ ~ Z^2®*P(-*'v/kT^) (1) 
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The source function for a two-level atom from Eq. 39, Appendix B, becomes 




2hV 

2 






k 


S 


where T represents the electron temperature and T the atom temperature, 
e a 

The expression for and as written above are for energies hv greater 

than the ground state ionization edge. For energies (hv) less than the ground 

state ionization edge, is zero. Likewise, for energies (hv) less than the 

excited state ionization edge, both and 02 ^ are zero. Also, is a function 

of frequency, because it includes the line shape. 

For the assumption that the excited states and electrons are in equilibrium 

at the electron temperature, one can write the expressions for Z^ 2 » ^2.f* 

Z in terms of T and T . The expression for Z,„ becomes 
2f e a 12 


’^=^2) rn (1)-| n u; g. 


n (2) 


(3) 


,(2)-x 


where x^2 the energy difference between the ground and excited state. Note that 
the bound electronic states are in equilibrium at the atom temperature {In ~ 


when the atom is in equilibrium. Since the excited state a^(2) is assumed to 


be in equilibrium with the electrons at T^» the Saha Eq. relates n^(2) to n^ as 


( 4 ) 


n 


n,(2) 


3/2 

. (2TTm kT ) . 

2 ^ .xpC-X/lcIe) . 

h"^ ®2 


( 4 ) 


where X 2 ionization energy for the excited state and gj^ is the ion degeneracy. 

If n (1) is in equilibrium with the electrons at the electron temperature, one has 
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( 5 ) 


n*(l) 

a 


2(27mi kT ) ' g. 


where ground state ionization energy. Eliminating n^ from Eqs. (4) 

it 

and (5), one can relate n (2) to n (1) and write 

a a .j, 

e 


n (2) n (l)'T' g, 

a _ a 

n (1) n (1) g. 


exp(-Xi2/^T^) . 


Thus, Z^2 written as 


^12 ' n (1) “'"‘kl T 

a a e 


where n (1) is given in terms of n and T by Eq. (5) 
am e e 

e 

approaches equilibrium Z approaches unity. 


Note that as the gas 


One can write as 


2 . 

r — — § r — „ 

'If n (1) 

a e 


"^e h h^ e xp(Xi/kTe) 
”a^^^ 8i 2(2TTm 


where the free states and the bound states are in equilibrium at the electron 
temperature when the gas is in equilibrium. If one eliminates n^ by using Eq. 


(5) , he obtains 


n*(l) 
a T 

'if “n^(D 
a 


where n (1)^ is defined in terms of the local electron number density and the 
e 

electron temperature by Eq. (5). Note that approaches unity as the gas goes 


to equilibrium. 



The parameter is always unity, because the excited states and electrons 
are in equilibrium at the electron temperature. 

As a consequence of Eqs. (3) through (9), one can rewrite Eq. (1) for the 


absorption coefficient as 

iia(l) 

k = n (l)a,,{l - — 7T\^ exp(-hv/kT )} + n (2)a_^{l - exp(-hv/kT )} 
V a If n V.1; e a 2f e 


+ n „( l ) a , -,{1 “ 


n*(l) V TT T '> 


iaV-/oi2X- - ^ (1) 

a 


exp[- 


kT T 
a e 


] exp(-hv/kT )} 
a 


( 10 ) 


Equation (2) for the source function becomes 


O _2h^ 
2 


[n (1) a + n (2)a 1 exp(-hv/kT ) 

a T ^ ^ 


V 


"a^^^T "^12 ®^P^-^l 2 ^^a " T^)/(kT^T^)] exp(-hv/kT^) 



( 11 ) 


As noted above, one must be careful in the application of these equations be- 
cause of the spectral behavior of and Note that as the gas approaches 

equilibrium, becomes the Planck function, and the absorption coefficient be- 
comes the product of the stimulated emission factor (1 - exp(-hv/kT)) and a 
term dependent upon the number densities of the bound states . The bound state 
population becomes Boltzmann distributed at equilibrium. 

Equations (10) and (11) show that the radiation source function and ab- 
sorption coefficient differ from their equilibrium values for the H 2 ”He loniza- 

(3 4) 

tion model. Tlwari and Szema * used equilibrium values for the source func- 
tion and absorption coefficient. They evaluated k^ and at the electron 

temperature. As discussed below, this assumption for k^ and leads to in- 
creased radiative heating for the outer planet atmospheric entries as compared 
to a total equilibrium solution. 
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FINITE RATE IONIZATION 


The finite rate ionization model for hydrogen-helium gas mixtures has been 
used in investigations of the radiative heating of vehicles entering the atmo- 
sphere of Jupiter and Saturn. Table 2 lists a few of these investigations and 
contains the ambient conditions, the probe nose radii, the shock wave stand-off 
distances, and the radiation fluxes at specific points in the entry trajectories 
for three cases. The solutions shown in Table 2 are for roughly the same entry 
conditions. Table 2 gives the comparison of the radiative flux for an equilib- 
rium solution and a finite ionization rate solution attained from the same 

( 2 ) 

numerical code. The solutions obtained with the Leibowltz code show that 

finite rate ionization produces a decrease in the radiative heating rate, 

whereas those obtained with the Tiwari and Szema^^’^’^^ code show that it 

produces an increase in the radiative heating rate relative to equilibrium 

ionization solutions. Note that in one of the cases, the equilibrium solution 

( 10 ) 

of Moss was used as a comparison case. 

The solutions of Moss^^^^ were obtained from an equilibrium viscous shock 
layer numerical code, which includes a detailed description of the radiative 
transport, equilibrium chemistry, and transport properties. The analysis al- 
lowed for ablation injection at the heat shield and either laminar or turbulent 
flow. The set of equations describing the flow were obtained from the steady- 
state Navier-Stokes equations by retaining terms up to second order in the in- 
verse square root of the Reynolds number. This numerical code was used to pro- 
duce "benchmark" solutions for the design of the Galileo Probe heat shield. 

(3 4) 

Tiwari and Szema * developed a viscous shock layer numerical code slmi- 

( 10 ) 

lar to that of Moss to solve the steady-state shock layer equations. Their 
analysis is restricted to nonblowing laminar flow fields about hyperboloid bodies. 
Equilibrium and finite rate chemistry and thermal equilibrium or non-equilibrium 
(only a two temperature gas - an electron temperature and a heavy particle temperature) 
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Table 2. 


Predictions of Equilibrium and Finite Rate Radiative Fluxes 



Conditions at infinity 
in front of shock 

Probe 

Nose 

Radius 

Shock 

Stand-off 

distance 

radi 

i 

Lation 

Elux 


km/sec 

T 

00 

K 

Poa , 
Kg/m^ 

% 

cm 

6 

cm 

« 

9eq 

Mw/m^ 

t 

Qnon-eq 

Mw/m^ 

Saturn Traj 
T = 51.75 sec 
atm = 89/11 

29.3 

83.9 

.000579 

31.1 

2.54 

40.0^ 

71. 7^^ 

(3) 

Saturn Entry 
Z = 115 Km 
atm = 85/15 

39.1 

145.0 

.000465 

23 

2.02 

422.6*^ 

1 

464.3^ 

Saturn Entry 
T = 18 sec 
atm = 73/27 

15 degree 

entry 

25 

i 


37. 0*^ 

20. 0^^ 


a. Moss, reference 10 

b. Te = Ta for non-equilibrium solution. Zoby and Moss, ref. 7 (using Tiwari and 
Szema Code, ref. 3) 

c. Tiwari and Szema, ref. 4, pp. 161 

d. Leibowitz and Kuo, ref. 2, Fig. 7 






















options are included in the code. The equilibrium chemistry uses an approximate 

equation of state. The finite rate chemistry uses the rate constants and ioniza- 

( 1 ) 

tion model of Leibowitz . A 58-step spectral model for the absorption coef- 
ficient of the hydrogen-helium gas mixture was used to calculate the radiation 
transport. 

( 2 ) 

Leibowitz and Kuo developed a simplified rapid numerical method to cal- 
culate a first order approximate solution for the flow field. They used cor- 
relations for the shock wave shape and stand-off distance and solved for the 
flow along stream tubes. Many details were omitted. The radiation calculations 
were uncoupled from the flow field calculations; however, they were correlated 
to the local thermodynamic state of the gas. A reaction similarity parameter, 

JJ, which represents the ratio of the fluid mechanic residence time to reaction 
time, was used to correlate the electron distribution in the shock layer. When 
was small 1), the ionization rate was finite. When Q was large 20), 
the ionization was in equilibritim. 

Leibowitz and Kuo Solution 

( 2 ) 

Figure 3 shows a shock layer solution obtained by Leibowitz and Kuo for 
the electron temperature and the electron concentration for two values of at 
a shock layer position slightly displaced from the stagnation point . For the 
equilibrium solution, the electron temperature is approximately constant across 
the shock layer. In the finite rate case, T^ is high just behind the shock and 
decreases to its equilibrium value near the body. For equilibrium flow, the 
electron concentration is roughly constant across the shock layer, and for finite 
rate ionization, it is very low near the shock and increases to its equilibrium 
flow value near the body. In the Leibowitz finite rate ionization model, the 
excited state population is related to the electron concentration. In other 
words, it would be small just behind the shock and increase to its equilibrium 
value near the body just as the electron concentration does . 
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ELECTRON CONCENTRATION 

mole/cm^xiO”® Tg, KxiO* 



Figure 3. Electron Temperature and Electron Concentration as a Function 
of Nondlmensional Position across the Shock Layer, = 20 
represents Equilibrium Flow, = 0.95 represents Nonequilibrium 
Flow, 6q = 1.7 cm. The Data is for 6 = 0.20 Radians, where 
0 = 0 at the Stagnation Point. (From Ref. 2) 
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Tiwari and Szema Solution 


(3 4) 

Tiwari and Szema ’ used a detailed numerical analysis to calculate 
the flow field and to predict the radiation heating of outer planet atmo- 
spheric entry vehicles. They calculated the shock layer composition by using 
the finite rate ionization model developed by Leibowitz. Figures 4 and 5 show 
a typical solution for the electron and atom temperature and the composition 
at the stagnation point for both equilibrium ionization and finite ionization 
rate flow. The solution has the same general trends as that of Leibowitz and 
Kuo discussed previously. 

Figure 4 shows that T and T for nonequilibrium flow are both greater than 
the equilibrium temperature just behind the shock, whereas, near the body the 
temperature is essentially the same for equilibrium and non-equilibrium flow* 
Thermal equilibrium occurs at about r| = 0*7 and the equilibrium and non-equi- 
librium solutions yield about the same temperature for values of r| less than 
0 * 2 . 

Figure 5 shows that the electron density increases very slowly behind the 
shock for finite rate ionization (nonequilibrltim) as compared to the equilibrium 
solution. This is an important difference between the two cases when radiative 
emission is considered, because the electron population is very important in 
the radiation transfer. Also, for finite rate ionization the atomic hydrogen 
population is much higher just behind the shock than the equilibritim hydrogen 
population. 

Figure 6 shows typical values of the shock layer optical thickness for the 
solutions shown in Table 2 as a function of photon energy hv. In general, the 
shock layer is optically thin in the spectral region from 1 to 9 ev. It is 
optically thick at the Lyman line at about 10 ev and for the Lyman continuum, 
at energies greater than 13.6 ev. The shock layer also becomes optically thick 
for free-free continuum radiation at very low values of photon energy. The 
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TEMPERATURE, KXIO 




MOLE FRACTION 



Figure 5. Species Mole Fraction Variation as a Function 

of Nondimensional Position across the Stagnation 
Shock Layer. (From Ref. 4) 
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OPTICAL THICKNESS 



Figure 6. Optical Thickness of the Stagnation Shock Layer as a Function 
of Radiation Energy, 
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lower part of Figure 6 shows the spectral position of the maximum of the Planck 
function as a function of temperature. For the temperatures in the shock layer 
the maximum occurs in the spectral region where the hydrogen shock layer is 
optically thin. Thus, when the radiation source function is the Planck function, 
most of the radiation emitted just behind the shock can easily reach the body. 

Figure 7 shows the equilibrium and finite rate ionization radiative heat 
flux toward the body as a function of nondimens ional position in the shock 
layer for the second case listed in Table 2. For the equilibrium ionization 
case the majority of the radiative emission occurs relatively close to the 
body at about q = 0.2. In the finite ionization rate calculations the most 
Intense radiative emission occurs near the shock wave at about q = 0.9. This 
difference can be explained by considering the radiation source function and 
the optical thickness in the shock layer. 

Radiative Flux at the Body 

The radiative flux reaching the body involves an integral over the source 
function (see eqs. (40) and (41), Appendix B) . 

qR(0,v) = -2l[ J J ^ S^(t^)exp[t^/y] dt^ dy. ( 12 ) 

-1 0 

Using the tangent-slab approximation, Eq. (12) can be integrated over y to yield 

/* 

q^(0,v) = -27ry ^ S^(t^)E 2 (Vdt^ (13) 

0 

where E 2 (t^) is the Exponential Integral of order 2. The source function 
S^(t^) and the optical depth are functions of the electron temperature 

and the electron and excited state populations as shown by Eqs (10) and 

(11) . When is small E 2 (t^ Is approximately unity and the radiative flux can 

be written as 
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i 


q^(0,v) ■ (14) 

0 

This relation approximates the radiative flux at the body in the spectral region 
from 2 to 9 ev, where the shock layers are optically thin. Equation (14) shows 
that the radiative fliox reaching the body is the integral of the source function 
over the optical radiating volume. It can be written in terms of the shock layer 
stand-off distance, 6^ as 

/* ^ 

q^(o.v) = -2Try s^(y)k^(y)dy. (15) 

0 

Thus, the radiative flux at the body is directly related to the radiation source 
function, the radiative absorption coefficient, and the shock wave stand-off 
distance when the shock layer is optically thin. 

Using Eq. (11) one can rewrite Eq. (15) as 

qj^(0,v) = -2 tt/* + n^(2)02^]exp(-hv/kTg) 

•/ Q i- T 

0 e 


+ n*(l) 




- V 

kT T 
a e 


]exp(-hv/kT^)J dy. 


(16) 


where it is assumed that the excited states and electrons are in equilibrium at 

T^. Care must be taken in considering Eq. (16) because and 

functions of the radiation energy. In the optically thin spectral region of the 

shock layer (approximately 2 to 9 ev) is zero. Also, one expects 

small in this spectral region because it is mainly in the line wing. Therefore, 

one can write the radiative flux as 

2 

qR^O*^) = -2 it f 

c 

0 
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Thus, the radiative flux is strongly dependent on T and n (2). The value of 

6 cL 

for hydrogen from page 265 of Reference 11 is 

^2f ^ 10”^^(hV2/hv)^ cm^ (18) 

where h \>2 is 3.4 ev for hydrogen and hv is the photon energy in ev units. 

Tiwari and Szema assumed the source function to be the Planck function in 

their analysis. In addition, the absorption coefficient was determined using 

equilibrium populations for the atomic electronic states. In their computer 

code the absorption coefficient was calculated from the number densities for 

atomic hydrogen, ionized hydrogen, and the electrons at the local temperature. 

( 12 ) 

They used a 58-step spectral model developed by Sutton , which automatically 
assumed that the atomic electronic states were Boltzmann populated at the given 
temperature. For the thermal nonequilibrium, finite ionization rate cases they 
used the electron temperature to evaluate the source function and absorption 
coefficient. These assumptions are equivalent to writing the absorption coef- 
ficient from Eq. (10) as 

= n^(2)o^^[l - exp(-hv/kTp] (19) 

and the source function from Eq. (11) as 


S 


V 



-hv/kT„ 

e e 

-hv/kT ’ 

1 - e ® 


( 20 ) 


where Eqs. (19) and (20) only apply to the optically thin region of the shock 
layer. Further, Eq. (17) for the radiative flux at the body can be rewritten 
using Eqs. (19) and (20) as 




(0 


,v) = 


-hv/kT 

® dy . 


( 21 ) 
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Equation (21) points out the strong dependence of the radiative flux on the ex- 
cited state population and the electron temperature. 

In the finite rate ionization solutions of Tlwari and Szema, is high 
behind the shock as can be seen from Fig. 4. This increases the exponential 
term of Eq. (21) . High values of T also increase n (2) , because they assumed 
that the electronic state population is Boltzmann distributed at T^. Thus, 
the finite rate ionization solutions produce a large radiative flux. 

In the equilibrium case the temperature is smaller just behind the shock, 
as shown in Fig. 4. This reduces the value of the exponential term in Eq. (21). 
In addition, the number of hydrogen atoms is much lower than it is for the finite 
ionization rate case (see Fig. 5) . Both of these effects reduce the excited 
state population. Consequently, the contribution of the gas near the shock to 
the radiative flux is reduced, compared to the finite ionization case. Near the 
body Tg and n^(2) are about the same for equilibrium and finite rate ionization; 
therefore, the radiative flux from this part of the shock layer will be about 
the same for both assumptions. Thus, the overall magnitude of the equilibrium 
ionization radiative flux is less than the finite rate ionization radiative 
flux . 

The discussion to this point has been limited to the optically thin ap- 
proximation and the Balmer region of the spectrum. However, the conclusions 
are generally valid. As can be seen from Fig. 6 the ground state continuum 
radiation is optically thick so only gas near the body contributes to radiative 
heating of the body in that spectral region. The atomic line transitions are 
also optically thick near the line center, which means that the spectral regions 
near the line centers will contribute to radiative heating of the spacecraft 
only if the radiation is emitted close to the body. The free-free continuum 
transitions occur mainly in the spectral range below 2 ev. The free-free con- 
tinuum radiation can be important; however, at the temperatures of interest its 
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contribution to the flux will be small, because of its small source function. 
Consequently, the radiative heating comes mainly from the Balmer region of the 
spectrum where the shock layer is optically thin and the above approximations 
hold. 

General Discussion 

(3 A) 

The physical problem Tiwari and Szema ’ solved is one in which the 
hydrogen atoms have a Boltzmann population distribution of their electronic 
states at the local electron temperature. The ionization rate was allowed to 
be finite by using the rate equations developed by Leibowitz^^^ . They in- 
terpreted the excitation rates in the ionization model as given in Table 1, re- 
actions 3, 4, 5 and 6, as ionization rates. Also, they did not assume the ex- 
cited states and electrons to be in equilibrium at the electron temperature 
as Leibowitz did. Leibowitz used the excitation rate as the ionization rate 
limiting step in a two step ionization model. Once the electrons were excited 
the number of electrons that became free was determined by the Saha equilibriinn 
condition at the local electron temperature. 

In the Leibowitz ionization model the excited state is populated at a finite 
rate. Thus, its population will lag its local equilibrium value. This implies 
that the excited state is not populated at its Boltzmann distribution value at 
the local temperature. Also, the excited state population and the electron 
population are assumed to be in equilibrium with each other at the electron 
temperature. Therefore, the excitation rate is not equal to the ionization 
rate. Both of these effects influence the radiation transfer, because the ex- 
cited state population distribution influences both the absorption coefficient 
(Eq. 10)) and the radiation source function (Eq. (11)). In addition, Tix\rari 
\ and Szema used the Planck function for the source function throughout the shock 

layer. As is pointed out above, the source function is equal to the Planck 
function only when the gas is in equilibrium. 
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Table 3 illustrates the difference in the population of the first excited 
state of hydrogen for the assumption of a Boltzmann population distribution at 
the electron temperature and for the assumption of the excited state population 
being in equilibrium with the electron population at the electron temperature. 

The data for T^, n^, and n^ just behind the shock wave was obtained from the 
solution of Tiwari and Szema shown in Figures 4 and 5. Note that for the Boltz- 
mann distribution the excited state population is much greater than it is for 
the assiimption that the excited states and electrons exist in equilibrium to- 
gether at the local electron temperature, 

(3 4) 

The radiation absorption coefficient subroutine used by Tiwari and Szema * 
automatically used the Boltzmann distribution for the excited states population. 
This greatly increased the absorption coefficient compared to its value for the 
assumption that the excited state population is in equilibrium with the electron 
population at the electron temperature. This effect together with the high elec- 
tron temperature just behind the shock compared to the temperature for equilibrium 
solutions (resulting in a larger value for the radiation source function) in- 
creased the radiative emission and is the reason that Tiwari and Szema calculate 


increased radiation in the finite rate ionization case. 


TABLE 3. Population of the first excited state of Atomic Hydrogen assuming: 
1) Boltzmann distribution at Tg and 2) excited state in equilibrium 
with electrons at Tg. (Tiwari and Szema solution) 


n 

T 

e 




n„(2)* 

** 

nH(2) 



ly 

cm 

1/cm^ 

l/cm^ 

1/cm^ 

1.000 

20,000 

8.39 

X 

10^’ 

4.66 X 10^^ 

7.93 X 10^^ 

9.14 X 10^ 

.950 

19,500 

8.49 

X 


9.99 X 10^^ 

7.67 X 10^^ 

4.59 X 10^ 

.900 

19,100 

8.22 

X 


3.64 X 10^^ 

6.57 X 10^^ 

6.57 X 10^ 

.875 

18,700 

8.26 

X 

10^' 

6.6 X 10^^ 

5.79 X 10^^ 

2.33 X 10' 

.850 

18,400 

8.16 

X 


8.62 X 10^^ 

5.17 X 10^^ 

4.21 X 10‘ 

.825 

18,000 

8.19 

X 


1.34 X 10^^ 

4.51 X 10^^ 

1.11 X 10' 

.800 

17,000 

8.20 

X 

10^’ 

1.53 X 10^^ 

3.08 X 10^^ 

1.77 X 10‘ 


n^(2) Boltzmann Population at T^ 
n,,(2) In equilibrium with electrons at T 


SUMMARY AND CONCLUSIONS 


Currently there is disagreement as to whether the cold wall radiative 

heating for probes entering the atmospheres of the major planets, increases or 

decreases relative to its equilibrium value when finite rate ionization is con- 

( 2 ) 

sidered. Leibowitz and Kuo predicted that finite rate ionization in the hy- 
drogen gas just behind the shock wave will reduce the radiative heating; how- 

(3 4 ) 

ever, Tiwari and Szema * asserted that it will increase the radiation heating. 

The present study investigated the radiation modeling used in each of the 
above mentioned investigations. It is shown that the radiation transfer analy- 

(3 4 ) 

sis used by Tiwari and Szema * over predicts the radiation emission because 
they assumed a Boltzmann population distribution at the local electron tempera- 
ture (which is higher than the equilibrixam solution temperature just behind the 
shock) for the population of the electronic states of atomic hydrogen. This 
equilibrium population assumption changed the nonequilibrium radiative source 
function to the Planck function. In addition, it also influenced the atomic 
radiation absorption coefficient for both line and continuum radiation through 
the excited state population. Both of these effects increased the local radia- 
tion emission just behind the shock wave in the 2 to 9 ev region of the spectrum, 
where the shock layer is optically thin. This, in turn, increased the radiative 
heat transfer to the body as compared to equilibrium chemistry predictions. 

The more realistic and generally accepted assumption is that the excited 
states and electrons exist together in equilibrium at the local electron tempera- 
ture. This is the assumption used by Leibowitz^^^ in the development of the hy- 
drogen ionization model. In the finite rate ionization model, excitation to the 
first excited state is the rate limiting step. Once the electron is excited it 
is ionized according to equilibrium conditions at the local thermod 3 mamic state. 
This implies that the excited states are not populated to their Boltzmann values 
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at the local conditions. 


Tiwari and Szema interpreted the finite rate excitation reactions in the 
ionization model as ionization reactions. If they would have correctly used 
the finite rate ionization model in their analysis they would have calculated 
smaller excited state populations and a smaller source function in the region 
just behind the shock wave. These effects would have reduced their finite rate 
ionization radiative heating predictions to values below those predicted by as- 
suming equilibrium chemistry in the shock layer. 

In the near future the finite rate ionization model as developed by Lei- 
bowitz should be correctly used in a detailed flow field and radiative heat 
transfer analysis like that of Tiwari and Szema. The results of such an in- 
vestigation need to be available as "benchmark” solutions. Simplified solu- 
tions which make use of the optically thin shock layer assumption and are com- 
putationally fast need to be developed and compared to the "benchmark" solutions. 
Once they are developed, the numerically fast, finite ionization rate schemes 
need to be upgraded to include the effects of ablation of the wall. The ablation 
products should be allowed to have finite chemistry and non-Boltzmann populations 
of their electronic and/or vibrational states. These effects should be included 
in predictions of the radiation heat transfer to probes entering the atmospheres 
of the outer planets. The radiative heating predictions need to be as accurate 
as possible, because the heat shield mass on these missions are very critical. 

One does not want to over predict the size of the heat shield because of the 
weight penalty. Yet if one under predicts the size of the heat shield the probe 
may burn up before it completes its mission. 
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APPENDIX A: REVIEW OF NONEQUILIBRIUM SHOCK WAVE STRUCTURE 


This section is devoted to a survey of the conclusions and results of 

theoretical and experimental studies of the relaxation phenomena in shock waves 

and the shock wave precursors. Relaxation phenomena behind shock waves have 

(11 13-17") 

been discussed in reviews and monographs ’ \ This section contains a 

discussion of the research conducted on the relaxation and precursors of non- 
equilibrium shock waves. This section does not contain an exhaustive account 
of all the research conducted on shock wave structure, but rather presents a 
discussion of the background experiments and theory leading to the present 
state of knowledge. 

Figure A-1 shows a schematic of the nonequilibrium shock wave structure 
for an atomic gas. It is typical of shock waves advancing at Mach numbers 
from 10 to 40, into a gas at 300 K and 1 cm Hg pressure. Just behind the shock 
wave there is a collision dominated relaxation region, wherein the shocked gas 
relaxes to its thermal and chemical equilibrium state. In the equilibrium 
region the gas cools due to the emission of radiation. Some of the radiation 
emitted behind the shock wave propagates into the region in front of the shock 
wave, where it is absorbed by the cold gas to form the precursor. The degree 
of excitation and ionization in the precursor is coupled to the relaxation re- 
gion through the radiation transfer. 

Relaxation Region Shock Structure in Atomic Gases 
Strong shock waves and high radiative heating rates that occur during the 
entry into the atmospheres of the outer planets have created a need to improve 
our understanding of the effects of the radiation, and collision processes in 
the relaxation process behind the shock wave. The hot gas radiative energy can 
propagate through the viscous shock wave and be partially absorbed by the cold 
gas in front of it. This absorbed radiation is referred to as trapped radiation. 
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Figure A-1. Nonequilibrium Shock Wave Structure. 


The absorption of this energy heats, excites and ionizes the cold gas; there- 
fore, the shock wave propagates into an excited and ionized media. This pre- 
ionization and preexcitation can influence the approach to equilibrium in the 
relaxation region behind the shock and the heat transfer to the entry vehicle. 
The radiation which is not absorbed by the cold gas is lost, and is referred 
to as radiation cooling. 

(18) 

Heaslet and Baldwin'" have obtained solutions for shock waves propagating 

in a perfect grey gas in which all of the radiation is trapped* The nature of 

their solution depends on the amount of radiative heating compared to viscous 

dissipation* If radiative heating exceeds viscous heating, it is possible to 

have a shock with no discontinuities. When radiative heating is less than 

viscous heating, the shock has a discontinuity produced by viscous dissipation 

(19) 

which is embedded within a larger inviscid region. Cohen and Clarke and 
Chow^^^^ have shown that this is a valid model of the shock wave, if the pene- 
tration length. of the radiation is large compared to the extent of the viscous 
region. Solutions for shock waves with combined radiation and viscosity with 

no restriction on the size of the penetration length have been obtained by 

( 21 ) ( 22 ) 

Traugott . Sen and Guess also solved the combined problem, but employed 
the Rosseland approximation to calculate the radiation transport in the shock 
wave. 


Radiation cooling effects have been observed in shock tube experiments by 
(23) (24) 

Petschek, et. al. , Redkoboradyi and Fedulov , and Oettinger and Bershader 

(25) (26) (23) 

. McCheshney and Al -Attar and Petschek, et . al . , made calculations 

of radiation cooling assuming the shocked gas to be optically thin. Solutions 

(27) 

including self -absorption effects have been obtained by Pomerantz and 

f 28^ 

Yoshikawa and Chapman'' ' . 

Ionization Relaxation 


Considerable Information is available in the literature on the relaxation 
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of atomic gases behind strong shock waves. Most of the analysis has been done 
for shock tube ambient conditions of about 0.001 atmosphere pressure, 300 K and 

Mach numbers from 10 to 40. Experimental results have been obtained using argon 

(25,29-35) (36,37) ^ ^ (1,38,39) 

, xenon and hydrogen as the test gases. Numerical 

results for the ionization relaxation are available for shocks in argon^^*^^~^^^ , 

(47,48) ^ , (49-53) ^ ^ (54,55) 

xenon , hydrogen and helium 

Ionization proceeds in a very complicated way behind shock waves in atomic 

(29) 

gases. Petschek and Byron showed that more than one ionization process was 
necessary to reach equilibrium ionization. If the electron concentration is 
sufficiently high, electron-atom collisions are the most probable ionization 
process. Petschek and Byron measured the electron-atom ionization rate. Their 
ionization rate was consistent with a two-step reaction scheme consisting of 
excitation followed by ionization of the excited state. 

* 

e + A A + e, 

* + 

e + A -*■ A + e + e. 

Furthermore, their measurements showed that the electrons were not in thermal 
equilibrium with the atoms during the ionization process. This occurred be- 
cause the electrons gained energy by collisions with heavy particles less ef- 
ficiently than they lost energy by creating ions. Thus, ionization by electron- 
atom collisions resulted in a net loss of energy in the electron gas, which then 
became colder than the atom gas. Thus, the gas behind the shock wave can be 
assumed to exist as two separate gases: the electron gas and the heavy particle 
(atom) gas. The electrons establish equilibrium among themselves at the electron 
temperature while the atoms and ions establish equilibrium among themselves at 
the atom temperature. Consequently, thermal nonequilibrium will exist in the 
relaxation region. 

During the early stages of ionization the electron concentration is too 
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small to support ionization by electron-atom collisions. Petschek and Byron 
speculated that reactions involving atom-atom collisions, radiation, and con- 
taminants contributed to the ionization in this region. Weymann^^^^ theorized 
that the most probable atom— atom reaction would be a two-step reaction similar 
to the two-step electron-atom reaction. 


A + A ^ A + A, 

* + 

A + A-»-A +e + A. 


This reaction was experimentally verified by Harwell and Jahn^^^^ . Impurity 
levels had to be reduced to a few parts per million before their effects be- 

/ CQ\ 

came negligible. The experiment was further refined by Kelly who confirmed 

and improved the results of Harwell and Jahn. 

The effects of contaminant reactions have been investigated in xenon shock 

(59) 

waves by Hacker and Bloomberg . They found a very complicated set of re- 
actions involving contaminant species and xenon atoms in various stages of ex- 
citation, and radiation. Contaminant reactions in argon were theoretically 
discussed by Morgan and Morrison^^^^ . 

Radiation can be important in the ionization process too. Bibermann and 
Yakubov^^^^ have shown that line radiation from the region following the region 
of nonequilibrium ionization excites argon atoms in the nonequilibrium ioniza- 
tion region. The excited atoms are easily ionized and contribute additional 
electrons, which reduces the time to reach equilibrium. Radiation can affect 
the ionization by creating electrons ahead of the shock front. Such precursor 

electrons have been observed in a re-entry by their effect on radar cross 
. (61-63) 


section 


and in shock tubes is described below. 


(49-52) 

In a series of papers, Skalafuris investigated the structure of a 

shock wave in atomic hydrogen. His analysis was done for conditions that exist in 
stellar atmospheres rather than shock tube, or planetary atmospheric entry. He 
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considered shock waves traveling at 30 to 70 km/sec into un-ionized atomic hy- 

Q —10 

drogen at a temperature of 5000 K, for a range of densities from 10 to 10 

3 -5 -4 

gm/cm , which corresponds to a pressure range of 4 x 10 to 4 x 10 atm. 

The results show the same trends as those in argon, in that the Lyman radiation 
is important in forming the precursor, the precursor raises the gas temperature 
behind the shock, the collisional relaxation region behind the shock is optical- 
ly thin to Lyman and Balmer radiation and that the electron temperature is less 
than the heavy particle temperature in the relaxation region. 

A general model for nonequilibrium shock wave structure in atomic gases, 
including both collisional and radiative processes, has evolved. It consists 
of a radiation induced precursor, and embedded viscous shock, an inner col- 
lisional relaxation region and an outer radiation cooling region as shown in 
Figure A-1. 

The relaxation mechanism to ionizational equilibrium occurs in two stages. 

The first is due to atom-atom collisions and the second is due to electron-atom 
collisions. Thus, one expects a region just behind the shock where the electronic 
state excitation and the ionization are produced by atom-atom collisions. When 
the number of free electrons becomes significant the electron-atom collisions 
rapidly drive the excitation and ionization to equilibrium. Throughout most of 
the relaxation region the heavy particles are at a higher temperature than the 
electrons because the electrons lose energy through ionization collisions faster 
than they gain energy by elastic collisions with the atoms. The nonequilibrium 
gas usually reaches thermal equilibrium just prior to attaining ionizational 
equilibrium. 

A model has evolved for the excitation and subsequent ionization of atoms. 
Most atomic gases have a large energy difference between the ground electronic 
state and the first excited state compared to the energy difference between the 
first excited state and ionization as shown schematically in Fig. A-2 for an 
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atom with two electronic states. The electrons and excited electronic states are 

very close to each other in energy; consequently, they are assumed to exist 

together in equilibrium at the electron temperature. In an experimental study, 

(9) 

Leibowitz, et. al. found that most of the excited atomic electronic states 
were in equilibrium with the free electrons and were at the free electron tempera- 
ture. 

It has been shown experimentally that the rate of ionization is controlled 
by the rate of excitation between the ground state and the first excited state. 

Once the atom is excited to the first excited state It will rapidly Ionize. 

Thus, the collisional ionization occurs as a two-step process as shown in Fig. 

A-2. In atomic gas shock wave structure studies it is acceptable to assume 

that the excited state population is in equilibrium with the electron population 

at the electron temperature and that the ionization rate is controlled by the 

excitation rate between the ground and first excited state. 

Nonequilibrium Shock Wave Structure in Hydrogen-Helium 

( 38 ^ 

Belozerov and Measures'' theoretically and experimentally investigated the 
initial ionization process in strong shock waves in hydrogen. Their analysis 
indicated that the electron temperature and atom temperature were essentially 
equal throughout the relaxation region. This appears to be in error due to the 
use of an incorrect value of the elastic atom-electron collision cross section^^^^ 
Leibowitz obtained results that agree with Belozerov and Measures when he in- 
creased the elastic electron-atom cross section to the value used by them. Thus, 
Tg should be less than T^ through most of the relaxation region as occurs in 
argon shock wave structure. Nakagawa and Wisler^^^^ measured the electron tem- 
perature in hydrogen shock waves and compared excited electronic state populations 
with predictions. They found that the higher excited states tended to equilibrate 
with the electrons at the electron temperature much faster than the lower excited 
states . 
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(Ij 9) 

Leibowitz ’ studied ionization rates behind shock waves in hydrogen- 
helium mixtures by measuring hydrogen line and continuum radiation. He developed 
a reaction scheme which included dissociation and a two-step excitation-ioniza- 
tion mechanism for hydrogen ionization by atom-atom and electron-atom collisions. 
He achieved good agreement between the numerical predictions and the experimental 
measurements. The excited states of hydrogen were assumed to be in equilibrium 
with the electrons at the electron temperature. The electron temperature was 
found to be significantly lower than the atom temperature throughout most of 
the relaxation region. 


Shock Wave Precursors 

The presence of electrons in front of hypersonic shocks has been observed 

(66-80) (81-102) 

and theoretically investigated by several persons . These pre- 

cursors are important in atmospheric entry because they influence the propaga- 
tion of electromagnetic radiation in the vicinity of the entry body. This 
affects identification of and/or communication with the vehicle. The precur- 
sor must be accounted for in interpreting the data a scientific probe gathers 
on atmospheric entry. Considerable effort has been spent determining the in- 
fluence of the precursor on heat transfer to atmospheric entry vehicles^^ 4,103)^ 
Precursor effects must be accounted for in the study of shock wave structure, 
since the precursor region is coupled to the relaxation region by the radiative 
transfer. The coupling of the precursor to the relaxation region has been stud- 


ied by several people 


(8,41,45,46,54,55,104,105) 


. In addition, astrophysicists 


must understand precursors because they occur in the shock structure of variable 

. (106) . ,, (107) . , 4 , (108) 

stars , interstellar gas , and certain solar events 

Precursors in Atomic Gases 

Throughout the literature controversy exists as to whether the precursor 
ionization observed in shock tubes is produced by electron diffusion from the 
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high pressure region behind the shock, or photoionization from the ultraviolet 
radiation emitted in the region behind the shock. Photoemission from the shock 
tube walls has also been considered as a possible source of the precursor. 

Many experiments have been undertaken to determine the mechanism for electron 
production. The earlier work leaned toward electron diffusion and photoemis- 
sion'' ; however, more recent experiments have conclusively shown that 

photo ionization of the gas ahead of the shock is the main mechanism for pro- 

^ . (72-76) 

ducing precursors 

Precursor photoionization may occur either as a one-step process wherein 
the electron is freed from the ground state of the cold gas in front of the 
shock by the absorption of ground state continuum radiation, or as a two-step 
process wherein the atom is first excited by line radiation and then ionized 


by excited state continuum radiation. These processes are shown schematically 

(85) (91) (92) (96) 

in Fig. A-2. Lagar'Kov and Yakubov , Murty , Vulliet , Dobbins 

and Nelson^^^^^ have investigated the importance of line radiation in the pre- 
cursor and the coupling of line radiation and excited state continuum radiation 
in the precursor. Precursors caused by photoionization gaseous impurities have 

( 74 ^ 

been investigated for trace amounts of hydrogen impurity in argon shock waves ’ 

90) ,(80) 

The electron temperature in the precursor has been measured 

( 66 ) 

The earliest precursor experiments were performed by Hollyer in a metal 
pressure driven shock tube. He investigated Mach 9 shocks in argon at pressures 
less than one mm Hg, using Langmuir probes . Far ahead of the shock his probes 
collected only electrons. Thus, he concluded that photoemission from the shock 
tube wall, due to radiation emitted from the shocked gas, was responsible for 
the precursor. 

Gloersen^^^^ experimented with xenon in a pressure driven pyrex shock tube 
at pressures from 0.75 to 4.0 mm Hg using probes wrapped around the outside of 
the tube. He observed two precursors, one traveling at about the speed of light 
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and one at the speed of the shock front. The fast precursor was attributed to 
photoemission from the shock tube walls, while the slow precursor was thought 
to be caused by collisional ionization, or photoionization of impurities at the 
shock front. Weymann, et. al. * ^ ^ used electrostatic and magnetic probes 

in a shock tube setup similar to that of Gloersen to investigate argon shock 
waves at Mach numbers from 8 to 12 and pressures from 2.5 to 10 mm Hg. They 
first concluded that the precursor was caused by diffusion of electrons from 
behind the shock wave because the electrostatic probe signals were negative 
and the magnetic probes indicated that electrons were moving ahead of the shock. 
They reasoned that the large concentration gradients just behind the sh6ck 
and high thermal velocities of the electrons were responsible for the electron 
diffusion. However, a later set of experiments at lower impurity levels showed 

(79) 

that photoionization was the cause of the precursor. McRae and Leadon also 
found impurities to be the main source of precursor ionization in their experi- 
ments in xenon. 

(72) 

Zivanovic investigated precursors with a ballistic range in which small 
copper spheres were fired up to 6,600 meters per second into air at pressures 
from 1 to 30 mm Hg. Electron density measurements were made utilizing a biased 
electrode plate, shielded by a screen. The bias established an electric field 
between the screen and the electrode plate. The field forced all the electrons 
which were photoemitted from the plate to return to the plate. Hence, only the 
electrons produced from diffusion and photoionization contributed to the cur- 
rent through the plate, which was measured and related to the electron density. 
For a series of experiments a second grid was placed in front of the screen and 
biased so that it would repel all the diffusion electrons. Thus, these experi- 
ments measured only the electrons generated by photoionization between the grid 
and the screen. The measured current was the same in the experiments which mea- 
sured only photoionization electrons and those which measured both diffusion and 
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photoionization electrons. Thus, it was concluded that the precursor was caused 
mainly by photoionization. 

Lederman and Wilson^^^^ measured argon precursor electron densities with 
a microwave resonant cavity and an electrostatic probe. The Mach numbers varied 
between 10.5 and 13.5 and the argon pressures were between 0.76 and 15 mm Hg. 

Some of the experiments were performed with a 'Venetian blind" in the shock tube. 
The "Venetian blind" allowed electrons to diffuse through it and blocked the 
radiation. Their results agreed with those of Zivanovic in that photo ioniza- 
tion, rather than electron diffusion, was the major cause of the precursor 
ionization. 

Precursors in Diatomic Gases 

The literature involving the influence of precursors in diatomic gases is 
usually Involved in an engineering application. In most cases it Involves the 
effect of the precursor on blunt body heating for planetary entry conditions. 

Precursors in diatomic gases have been observed ^ and theoretically 
(93,94,99 103) 

modeled » » > ^ Much of the radiation energy emitted in the shock layer 

relaxation region is lost; however, when the photon energy is greater than the 
molecular dissociation energy, it is strongly absorbed in the precursor region. 
This absorbed energy dissociates and ionizes the gas. It also increases the 
gas temperature and pressure. The change in the free stream flow properties 
may change the flow characteristics of the shock layer and the heat transfer 
to the body. 

(75) 

Mermagen observed precursors in front of Lexan cylinders, launched from 
a light gas gun at velocities from 3900 to 5500 m/sec into air at a pressure of 
50 mm Hg. The excited and ionized gases in shock layer around the model emitted 
Infrared radiation which pre-excited the cold gas in front of the shock. The 
shock layer was populated by excited and ionized air, including contaminants, 
as well as ablation products from the Lexan. Lin^^^^ presented evidence of a 
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"photo ionization halo" surrounding the first manned Mercury orbital flight during 
its re-entry. 

Radiative preheating of the cold gas by absorption of photons emitted from 
the hot shock layer was considered in most shock tube studies of shock wave 
structure by the late 1960 *s; however, it was usually neglected in the analysis 
of the blunt body shock layer problem. Yoshikawa^^^^^ considered the precursor 
influence on the blunt body shock layer for a gray gas and a linear relation- 
ship between the Planck emissive power and the gas enthalpy . Lasher and Wil- 
(103) 

son considered an uncoupled precursor, in which they used available blunt 

body shock layer solutions to obtain a first approximation of the radiative 

energy available to form the precursor. The precursor effects then preturbed 

the free stream conditions and changed the shock layer solutions. Thus, their 

solutions were valid when the precursor effects were small. Lasher and Wilson 

included the spectral variation of the radiation and radiation properties. 

(97) 

Liii investigated precursor effects on the heating of vehicles entering the 

earth's atmosphere in terms of the ratio of radiative to convective fux. The 

precursor influence was negligible tjhen the radiative to convective flux ratio 

was less than unity. 

(93) (94) 

Smith and Edwards calculated electron densities in the precursor 
of a shock preceding a blunt body for earth atmospheric entry. Smith used a 
simple microscopic model for air and included three-dimensional effects by 
modeling the body and shock layer as a point source. He calculated velocity, 
pressure, density and electron concentrations in the precursor region. Edwards 
used a more descriptive microscopic model for air and modeled the blunt body 
shock layer as a constant temperature disk. A constant fraction of the radia- 
tion absorbed in the precursor was assumed to yield the electron concentration 
along the stagnation streamline. The influence of the precursor on the blunt 
body heating for lunar returns was predicted to be negligible. 
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In the past 10 years there has been considerable interest in designing 

probes to enter the atmospheres of the outer planets, like Jupiter and Saturn. 

These atmospheres are thought to be composed of mainly hydrogen and helium. 

Thus, considerable effort has been put forth to investigate the heating of 

vehicles entering hydrogen-helium atmospheres. As a part of this effort Tiwari 
(3-4) 

and Szema have investigated the effect of precursors on the heating of 

vehicles traveling in hydrogen-helium atmospheres. They solved the radiation 
coupled precursor-shock layer problem taking into account many details such 
as viscous flow, spectral details of the radiative transfer, thermal non- 
equilibrium in the shock layer, multi-dimensional flow fields, finite chemical 
rates, etc. Their results generally show that the precursor effects should in- 
crease the heat transfer to the entry vehicle. 
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APPENDIX B; FORMULATION OF RADIATIVE INTENSITY 


This section of the Appendix contains a general formulation of radiation 
intensity. The medium through which the incident radiation field passes is as- 
sumed to be an atomic plasma. The plasma is assumed to consist of ground and 
excited state atoms, ions, and electrons. In formulating the radiation intensity, 
spontaneous and induced emission and absorption will be considered. This 
model makes it possible to study the coupling of line radiation (bound-bound 
processes) with continuum radiation (bound-free processes) . 

The fundamental quantities of atomic radiation are the Einstein coef- 
ficients^^l The Einstein coefficient, A^j^, for spontaneous emission is 
defined such that the probability per unit time for a spontaneous transition 
from upper state u to lower state % with the emission of radiation in the 
solid angle dQ is 

A - dn, (1/sec) (1) 

u5, 

The total number of spontaneous transitions from state u to per unit volume 
per unit time with the emission of radiation in dJ^ is 

n (u) A . dfi, (l/(cm^-sec)) (2) 

a ux 

where n (u) is the number density of the upper excited state. The energy 
a 

emitted per transition is equal to hv^^^ when the influence of the surrounding 
plasma is negligible; however, in real situations the transition can occur 
with the emission of almost any amount of energy. A certain value of emitted 
energy is determined by the emission line shape, (sec) which essentially 

represents the influence of the neighboring gas on the atom of interest. The 
total radiation intensity emitted in the frequency interval V to v + dv per unit 
time is 

hv n,(u) A „ ^ dQdv (ergs/(cm^-sec)) . (3) 

a uJl V 
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Spontaneous emission also occurs when a free electron recombines to a 
bound state. In the model, there are several possible free-bound transitions: 
free to first excited (f to 2) and free to ground (f to 1), etc. 

The probability that an electron-ion encounter occurs with the free elec- 
tron losing energy, hv = > by emitting radiation within the solid angle 

dJ^ is dfi. After the encounter, the electron-ion pair becomes an atom 
with the electron in state i. For the model, i can be any bound electronic 
level. In this notation Xjj Is the ionization potential of the £th electronic 
level of the atom and H is the translational energy of the free electron. 

Spontaneous emission cannot occur unless an electron and ion undergo 
an encounter. The number of collisions between electrons with energy in 
the range ri, n + dri and ions per unit time per unit volume, is 

A „ n n. f (ri) dridJ^ (l/(cm^-sec)) (^) 

fJl e 1 e 

where f^Cri) is the energy distribution function for the electrons. For a 
Maxwellian distribution at the electron temperature T^, 

f(ri)dn = 27 T(iTkT )~^^^(n)^^^exp(-n/kT )dn. (5) 

The probability of a spontaneous emission, with the creation of an atom, 
is equal to the probability of an emission times the number of electron-ion 
collisions. The total radiation intensity emitted in the frequency interval 
V to V + dv per unit volume per unit time in the solid angle dfi by all electron- 
ion spontaneous emissions is 

hv A-„ n n f (n) dridf^, (ergs/ (cm -sec)) (6) 

rS, e i e 

where dri = hdv. This process emits continuum radiation; whereas, the pro- 
cess of Eq. (3) emits line radiation. 

The Einstein coefficient for absorption is defined such that in a 
radiation field of specific intensity, I^, the probability per unit time for 
a transition from state 2. to u with the absorption of radiant energy in the 
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frequency range v to V + dv in the solid angle dfJ is 

B^u K \ (1/sec) (7) 

where is the absorption line shape. The total number of transitions from 
state to u per unit volume per unit time as a result of the absorption of 
radiation is 

na<« dfidv. (l/(cm^-sec)) (8) 

The total amount of radiation absorbed from the radiation field in the fre- 
quency interval V to v + dv per unit volume per unit time in the solid angle 
dQ is 

hv n^(Jl) B^^ (})^ dfidv, (ergs/(cm^-sec)) (9) 

This process represents a bound-bound process. The line shape, is used 
because the probability of the absorption of energy is related to the magnitude 
of the energy, through the line shape. 

The bound-»free process must also be considered. Equation (9) can be used 
to represent the bound-free process if the subscript u is changed to f. The 
continuum absorption can then be written as 

hv n^OO dfidv. (ergs/(cm^-sec)) (10) 

where % represents any bound electronic state. 

Emission of radiation may also be induced by the radiation field. The 
Einstein coefficient for induced emission is defined such that in a radiation 
field of specific intensity, I^> the probability per unit time for a transition 
from state u to as a result of induced emission in the frequency range v to 
V + dv and in the solid angle dfi is 

B , I dfidv. (11) 

ui V 
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The total number of transitions from state u to £ per unit volume per unit 
time as a result of induced emission is 

n^(u) dR, (l/(cm^-sec)) (12) 

The total amount of radiation emitted in the frequency interval v to v + dv 
in the solid angle dJ^ per unit volume per unit time by induced emission is 
then 

hv n„(u) B 41 I df^dv, (ergs/ (cm^-sec) ) (13) 

O’ UX/ V V 

where again the line shape, is introduced because it is the probability 

density function of the magnitude of the emitted energy. 

Equation (13) represents bound-bound induced emission. The induced emis- 
sion which results in continuum radiation must also be considered. Induced 
emission, when the initial states are free states, depends upon the number of 
collisions between the electrons and ions. Using Eq. (4) and Eq. (13), one can 
write an expression for the continuum induced emission. The radiation emitted 
in the frequency interval, V to v + dV, in the solid angle dS^ per unit volume 
per unit time by free-bound induced emission is 

hv B-. n n f (r|) I dndf2, (ergs/ (cm^-sec) ) (14) 

I X/ G 1 G V 

where dq = hdv. During the process, the electron loses energy, hv = + U, 

while the ion gains energy, Xjj^> recombination is at the £.th level of the 
atom, as the electron- ion pair create an atom in state £. 

The rate of change of the radiation intensity at a point within the solid 
angle, df2, in the frequency range, v to v + dv, in the direction s is the sum 
of the contributions (both negative and positive) given by Eqs. (3), (6), (9), 
(10), (13), and (14): 
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(15) 


dl 


dvdJi = {-[n (il) B + n (£) B <t> - n n, B f (n) 

ds a If £u V e i fz 


- V '•'v' '"’'v K 


+ n n A f (n)] hv} dvdfi. 
e 1 rio e 


The processes represented on the right-hand of Eq. (15) in the term multiplied 
by are respectively continuum absorption, Eq. (10), line absorption, Eq. 

(9), continuum induced emission, Eq. (14), and line induced emission, Eq. (13). 
The second term, which is multiplied by hv represents line spontaneous emission, 
from Eq. (3) and continuum spontaneous emission, from Eq. (6), respectively. 

Equation (15) represents the general form of the radiation transfer equa- 
tion. It can be rewritten as 


dl\) n n. B 

^ - {n(z) B [1 - f (n) 

n (u) 

+ n(Z) B„ (f... [1 - hvl 


a - iu % 
n-(a)A.o t 


n^(0 Ba (j) 
a Jlu V 


n 




a 




\f "e' 

a 


(16) 


hv. 


if 


I = 1,2,. . . L-1 
u = 2,3, . . . L 

where L is the total number of bound electronic states. In what follows the 
absorption and emission line shapes will be taken as the same, = <}>^. 

The relationships between the Einstein coefficients can be found by con- 
sidering the equilibrium case. If the gas is in equilibrium, = 0 and I^ 

is given by the Planck function 


2hv^ 

»V="T" 

C 

Considering only the line transitions 
one can write at equilibrium 


1 

exp(^) - 1 

, those terms involving ij)^ in Eq. 


(16) 


(17) 
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1 


( 18 ) 




n n (u)B . 

a In ^ a u£ 

^ ■ n,«TB7' 

a Jlu 


However, at equilibrium 


n*(u) 

a 


~ “P<- 


(19) 


where and g are the degeneracies of the states £ and u respectively, T 
£ u “ 

is the atom temperature, and X is the energy difference between the states. 

X/TJL 

To force Eq. (18) to agree with the Planck function one must have 


u£ 

^£u 


’£ 2hV 


B 




and 


u£ 

^£u 


u 


(20a) 


or that 


. 2hv^ '■n.(J)- 


£u 


rn (£)-| 

fcl 


hv 


exp(- 


£u\ 


kT 


^ _ rn*(«- 


B 


£u 


rn 

_ a 
-n (u)J„ 


hv 


exp(- 




kT, 


(20b) 


Then one can write 1 . as 


^v = 


2hv 


3 In^(u)/n^(£)] [n^(£) /n*(u) ]q; exp(- hVjj^^/kT^) 


c^ 1 - [n^(u)/n^(il)] tn^(£)/n^(u)lT exp(- hv^^kT^) 


(20c) 


which approaches the Planck function as the gas goes to equilibrium. This 
then defines the three Einstein coefficients, if one of them is known. 

When the continuum portion of the intensity is considered, one must define 
the Einstein coefficients as 


3, 2 


^f£^®£f ~ ®f£^®£f ~ 8jj^/(g£g^) 


( 21 ) 


where g^ is the degeneracy of the free electron states with energies between 


n and Ti + dn; 


(111) 


«f 


2(47t) (mg) 


3/2 


(2n)^^^ do. 


( 22 ) 
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The factor is the degeneracy of the Hth bound state, and g. is the degen- 

X/ X 

eracy of the ion. The ion is normally in its ground state. In Eq. (22), g^ 

3 

has units of 1/cm , so that it represents the degeneracy per unit volume. One 

^e^i 

must consider the relation — "i; — which appears in Eq. (16), where 

e 


f„(n) dn = 


2tt 


(irkTg) 


3/2 


(Tl)^^^exp(- 11 /kT ) dn 


(23) 


represents the fraction of electrons with energies between n and n + dn, where 
hv * Xj^ ^ . The relation can be written as 

exp(- n/kT^) 


e i 

n (^>) 
a 


2g,. <2TrmekT^) 


3/2 


using Eqs. (21), (22), and (23). At equilibrium the ratio of the population 

of the electrons and ions to that of bound state 5. is given by Saha equation 

,3/2 


(4) 


* * 
n n 

[-4-^1 


n*(Jl) T ® 

a e 


g. (2irm^kT^) exp(-Xj^ZkTg) 


(24) 


t 


Thus, the relation can be written as 


n n n*(j!,) 

exp(-hvZkT^) 

^ “i% 


(25) 


where the definition of 1) was used (hv = Xj^ • Using Eq. (25), one can 


write the continuum intensity as 

n n n*00 

[— § — y] exp(-hv/kT ) A /B 

"i T e f£ if 

a e X a 


I = 
V 


n rT ^;i) 

1 - “^TTT t 

n^U) n^n. ^ 


(26) 


exp(-hv/kTg) 


Now as the gas goes to equilibrium and using equation (21) for 

I -> B , the Planck function. Note that Eq. (26) is a function of the electron 
V V 
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temperature T^. 

Using the definitions given by Eqs. (20), (21), and (25), one can rewrite 
Eq. (16) as 

^ = - {n,(£)B^^[l - Z^^exp(-hv/ky ] + 


[1 - Z^^exp(-hv/kT)]}hvI^ 


+ {n^(£)Bj^^Z^^exp(-hv/kTg)hv + n^(S,)B^^4>^Z^^exp(-hV/kT)hv} 


2hv' 


(27) 


where 




n^Hj r n*(JQ 1 

■ -a(^> LnV J, 


e 1 T, 


(28) 


and 




n (u) rn (il)*l 
« a I a 

n ( 1 ) */ s 

a Ln (u)J 


a' ' T 


(29) 


and a is summed over all bound levels it » 1,2. . . L and u is summed over 

levels £ + 1, £ + 2 , . . . L. If the gas is in equilibrium, and Z^^ go to 

unity and dl^/ds goes to zero. The general form of the equation of radiative 

( 110 ) 


transfer is 


ds " % 


(30) 


where k^ is the absorption coefficient, and e^, is the emission function. From 


V 


Eq. (27), one sees that 


k^ = (n^(Jl)B [1 - Z exp(-hv/kTg)l 


+ ■ Z^^exp(-hv/kT)]}hv 


(31) 


and 


2hv 


T 4> Z exp(-hv/kT]hV (32) 


B-8 


In the limit of equilibrium, and for the lines agree with relations given 

by Jeff eries^^^^^ when the absorption and emission profiles are the same. The 
radiation transfer equation is simplified by converting the independent vari- 
able from s to T^, the optical depth. 


T 


V 



(33) 


The equation of transfer can then be written as 

dl 


y 


V 


dx 




(34) 


where is known as the source function, and y = ps, where y is 
the cosine of the angle between the general direction of radiation intensity 
propagation s and the coordinate y. The source function becomes 


hv hv 





(35) 


The radiative transfer can now be evaluated, if the line shape, (j>^y is known 
and the Einstein coefficients, and determined. and B^^ 

are functions of the particular gas of interest. Note that in the limit of 

equilibrium, reduces to the Planck function. Note also that the units of 

2 2 
are cm /erg and that the units of B^^ are cm /(erg-sec) because must 

2 

have units of cm /erg. 


as 


The absorption coefficient can be written in terms of the cross sections 

hv 




hv 


£f Jlf 


(36) 


where the cross sections are defined as 




(37) 
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and 


"Ju ■ '"’Vlu > 


(38) 


In terms of the cross sections, the source function becomes 


_ 

VT 

V 2 


hv 

kT 


hv 


hv 


n,(«a„(l - + njwojl - 


(39) 


The cross sections are functions of the gas and the various transitions. The 

general solution of the equation of transfer for the geometry shown in Fig. B-1 

( 110 ) 


is 


I (t ,y) = 

V V 


1 , 1 , 


(AO) 


0 "^v/^ r ^ 

(0,y)e + J S^it^)e — (u>0) 


where (tJ* ,y) is the intensity incident on the boundary at , and 

(0,y) is the incident intensity on the boundary at zero. The contribution 
to the radiative flux in the frequency range V to V + dv is 




(41) 
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